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Quantum mechanical calculations are performed of the emission and absorption profiles of the 
■ lithium 2s-2p resonance line under the influence of a helium perturbing gas. We use carefully 

' constructed potential energy surfaces and transition dipole moments to compute the emission and 

absorption coefficients at temperatures from 200 to 3000 K at wavelengths between 500 nm and 
1000 nm. Contributions from quasi-bound states are included. The resulting red and blue wing 
profiles are compared with previous theoretical calculations and with an experiment, carried out at 
a temperature of 670 K. 
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I. INTRODUCTION 



It has been argued that prominent features in the spectra of brown dwarfs and extra-solar planets may be attributed 
to the resonance lines of the alkali metal atoms, broadened by collisions with the ambient hydrogen molecules and 
helium atoms [1, 2, 3, 4, 5]. The interpretations of the measured profiles which yield information on the temperatures, 
^ , densities, albedos and composition of the atmospheres are based on models of the line-broadening. Recent studies of 
O ■ the potassium and sodium lines have employed classical and semi-classical scattering theories [5, 6, 7]. Jungen and 
Staemmler [8] have obtained detailed emission and absorption profiles (in arbitrary units) for the resonance line of 
lithium for a wide range of temperatures and Erdman et al. [9] have presented measurement of the absorption profile 
i-C ' of lithium in a gas of helium at a temperature of 2033 K together with the results of a quasistatic semiclassical theory. 

In this paper, we consider the pressure broadening of the 2p-2s resonance line of lithium arising from collisions 
with helium atoms. We carry out full quantitative quantum-mechanical calculations of the emission and absorption 
coefficients in the red and blue wings. We include contributions from quasi-bound states and we allow for the variation 
of the transition dipole moment with internuclear distance. We present results for temperatures between 200K and 
3000K and we compare with the emission coefficient that has been measured at 670K [7] and with calculations of 
■ Herman and Sando [10] at 670K. 
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For a system consisting of a mixture of lithium in a bath of helium atoms, the emission spectrum of the lithium 
2p-2s atomic line has both blue and red wings due to collisions with helium. If the gas densities are low enough 
that only binary collisions occur, the problem is reduced to the radiation of temporarily formed LiHe molecules. The 
broadened 2p-2s atomic emission line corresponds to transitions from the excited A 2 n and B 2 E electronic states of 
the LiHe molecule to the ground X 2 E state. The X 2 E and B 2 E states have no bound states and the A 2 n state has a 
shallow well which supports seven bound ro-vibrational levels. The level v = 6 is very close to the dissociation limit 
and it may not be found in practice. In any case it behaves as if it belonged to the continuum. We consider both 
bound- free and free- free transitions. 

The emission coefficient is defined as the number of photons with frequency between v and [y + dv) emitted per 
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unit time per unit volume per unit frequency interval, which for bound-free transitions is given by [10, 11, 12] 
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and for free-free transitions by [10, 11, 12] 
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where u> is the probability that the initial excited electronic state of LiHe is populated by the collision, which is | for 
the A 2 n state and | for the B 2 E state. In Eqs. (1) and (2), nu and njj e are the number densities of lithium and 
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helium atoms, respectively, c is the speed of light, and Qt = (2ir[j,kT/h 2 ) 3 > 2 is the translational partition function for 
relative motion, /i is the reduced mass of the collision pair, k is Boltzmann's constant and T is the temperature. T b f 
and r^f are defined by 

rb/ =EE^ J )K3f b ' J |^l^ J )| 2 , (3) 

V J 

and 

r" = £(2J+ 1) I dE a eM-E a /kT)\{g^ J \D\g^ J )\\ (4) 
J J 

where g(v, J) is the relative population of levels (v, J) of the A 2 n state and D{R) is the electronic transition dipolc 
moment which varies with the nuclear separation R. The subscript a denotes the initial state and b the final state. 
The functions g v ' J and g E ' J are, respectively, the bound and energy normalized free wave functions determined from 
the radial Schrodinger equation 

_____ + ¥ [E- V(R) 2tiR2 ]g(R) 0. (5) 

The bound energy E^' J and free energies E a and Ef, are measured from the dissociation limit of the excited states 
and the ground state respectively. The photon frequency v is determined by the relation hv — E a — E^ + hvo, where 
v n is the atomic line frequency. In writing Eqs. (3) and (4), we may replace the sum of the matrix elements in which 
J changes by ±1 by matrix elements in which J docs not change. 

The sum of Eqs. (1) and (2) is the total photon emission rate. In the limit of high densities, the bound levels of 
the A 2 IT state are in thermal equilibrium and g(v, J) is given by 

g(v, J) = (2J + 1) exp(- E%< J /kT). (6) 

In the limit of low densities, no bound states are populated and only free-free transitions contribute to the emission 
spectrum. 

For absorption, we consider free-bound and free- free transitions from the ground X 2 E electronic state to the excited 
A 2 II and B 2 E states. The absorption coefficient for free-bound transitions is given by [13, 14, 15] 

al~ b (b -f a) = ujULinHe^^-^^eMH" ~ vo)/kT] ■ T b f , (7) 
and for free-free transitions by [13, 14, 15] 

i f lGlT^V 1 7T r f 

al- f (b -> a) = um Ll n He — — — -3 exp[/i(i/ - vo)/kT] ■ . (8) 
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In Eqs. (1), (2), (7) and (8), the emission and absorption coefficients are given in terms of the number of photons 
in a frequency interval dv. The photon energies emitted or absorbed are obtained by multiplying k v or a v by hv. 
The coefficients per unit wavelength are obtained by multiplying by v 2 /c. 



III. POTENTIALS AND DIPOLE MOMENTS 



The potential energy surfaces were constructed to reflect recent calculations. For short and intermediate internuclear 
distances R, the energy points were taken from Ref. [16] for 3.0 < R < 8.0, Ref. [17] for 9.0 < R < 19.0 and Ref. [18] 
for R = 20.0 for the X 2 Y> state, from Ref. [18] for 2.0 < R < 2.75 and Ref. [19] for 3.0 < R < 20.0 for the A 2 H state, 
and from Ref. [20] for R = 2.0 and Ref. [18] for 3.0 < R < 20.0 for the B 2 Y, state. At large R, the potential energies 
vary as ~C 6 /R 6 . The theoretical values of C 6 of 22.51 for the A 2 E state [21], 28.27 for the A 2 II state [22] and 50.69 
for the £? 2 S state [22] were adopted. For each potential energy curve, the data in different segments were smoothly 
connected with a cubic spline procedure. The potential energy curves are shown in Fig. 1, along with the difference 
potentials V^2 n (_R) — Vx2_(i?) and Vb^s(R) — Vxz?,{R)- In Fig. 2, we present the energy separations as wavelengths 
A obtained with the classical assumption that the emission occurs through vertical transitions as the particles move 
along the potential energy curves. The vibrational energy levels of the A 2 ±l state are listed in Table I. Values of 
the difference of energy of the v — and 1 vibrational levels of the A 2 ±l state have been derived from spectroscopic 
experiments [23]. The measured values are 289.9(4) cm -1 and 290.0(3) cm" 1 , obtained by differencing two A 2 IT-D 2 A 
transitions. From Table I, we obtain 288.4 cm -1 . An earlier calculation by Jungcn and Staemmler [8] yielded 260 
cm -1 . 

The adopted dipole moment curves [24] are shown in Fig. 3. 



3 



IV. CALCULATIONS 



The radial equation (5) was solved using the Numerov method. The integrands of the free-free matrix elements 
oscillate at large internuclear distances and the integrals fail to converge. They can be transformed by the equation 

D(R) = {D(R)-D(oo)} + D(oo). (9) 

The matrix element of {D(R) — D(oo)} converges because the operator tends to zero at large R. The matrix element 
of the constant D(oo) may be written in the form[10] 

D(oo)( 9 ^ J \g^ J ) = D(^)(g^ J \AV\g^-')/(E a - E b ), (10) 

where AV = V a (R) — Vb(R) is the difference potential. The right hand side of (9) is convergent since AV — > for 
R — > oo. 

The integration over energy in Eq. (4) was carried out using the Gauss-Laguerre method with 100 points. We also 
ran the code with 200 points and found no difference in the results. We included values of J up to 50 for T = 200 K 
and up to 200 for T = 3000 K. 

Shape resonances occur in scattering by the A 2 II potential. They were determined by calculations at high energy 
resolution. Table II is a list of the values of J and the locations and widths 7^' J of the quasi-bound levels. We used 
Simpson's rule with energies selected at closely spaced intervals to evaluate the resonance contributions to the energy 
integrals. 



V. RESULTS AND DISCUSSIONS 



The bound- free emission coefficient (1) is a weighted sum of the emission rates of the individual ro- vibrational levels 
of the A 2 II state. The spontaneous transition probabilities are given by 

In Fig. 4 we plot W v in atomic units as a function of wavelength for rotational states J = 0, 5, 10 and 18 of the 
vibrational levels v = to v = 5. The corresponding lifetimes in seconds are listed in Table III. For high vibrational 
levels, they are approaching the limit of the radiative lifetime of the resonance state of the lithium atom which is 
27.2 ns [25]. 

In Fig. 5 we present the contributions to the emission coefficient k v in units of cm _3 s _1 Hz _1 of the A 2 II state from 
the bound-free transitions, assuming that the levels are populated in thermal equilibrium at temperatures ranging 
from 200 K to 3000 K. The contributions from the resonances were multiplied by 7"' J /(7r' J + lT J )> wnere 7r' J 
and ji' J are the resonance and radiative decay widths[ll], but the effect is negligible. At no temperature and at 
no frequency is the quasi-bound contribution more than a few percent of the total. The shape of k v results from a 
superposition of the spectra in Fig. 4 weighted by the populations of the ro-vibrational levels. The sharp increase 
near the resonance wavelength at 670.8 nm is the contribution from low-lying rotational states of vibrational level 
v = 5, illustrated in Fig. 6. A peak in the spectrum at around 870 nm is apparent at low temperatures but it tends 
to smooth out at high temperatures. At the higher temperatures, there is a plateau region associated with emission 
from the attractive region of the potential. The emission decreases rapidly at long wavelengths because there are no 
values of the internuclear distance R that correspond to wavelengths beyond 1000 nm, as shown by Fig. 2. 

The calculated free-free photon emission rates k v in units of cm _3 s _1 Hz _1 produced in transitions from the A 2 II 
state to the X 2 E state are presented in Fig. 7. Weak quantum oscillations are found that correspond to transitions 
from the attractive region of the A 2 II potential. There occur rapid increases in k v at wavelengths close to the atomic 
resonance wavelength. At some point as we approach the line center the binary approximation that we have adopted 
breaks down. Our calculations are valid only in the wings of the line. Studies of the profile at line center have been 
reviewed by Lewis [26]. 

The blue wings arise from the B 2 S-X 2 S transitions. Only free- free transitions occur and the result is a smooth 
profile decreasing rapidly with decreasing wavelength except for a weak satellite feature near 536 nm that is apparent 
at high temperatures. The behavior is expected from a consideration of Fig. 2 which shows that internuclear distances 
in the B 2 S state that correspond to short wavelengths are inaccessible. The satellite feature is produced by emission 
from internuclear distances at which the two potential energy curves are parallel [5] . Previous work has located it at 
516 nm [20] or 506 nm [7]. There is an indication of the satellite near 530 nm in an experimental spectrum obtained 
by Lalos and Hammond [27]. 
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Semi-classical calculations by Jungen and Stacmmlcr [8] over a range of temperatures above 200 K have similar 
profile shapes to those we obtain with quantum-mechanical calculations. Measurements of the line profile at 670 K 
at high pressures have been carried out by Scheps et al. [7] and calculations at high and low pressures in the red 
wing have been made by Herman and Sando [10]. In Fig. 8 we give a comparison with our results. There is broad 
agreement between the theoretical low pressure profiles which arise from the free-free transitions. The differences in 
detail can be ascribed to the adopted interaction potentials. The theoretical calculations of the high pressure profiles 
have the same shape. They both have a plateau at intermediate wavelengths, increase sharply near the line center 
and fall off rapidly far in the red wings. They differ quantitatively in the plateau region. Agreement between the 
present theoretical red and blue wings and experiment is close and it indicates that the binary approximation is valid 
to within 20 nm of line center. There may be a discrepancy at wavelengths longer than 950 nm where the emission 
coefficient is becoming very small. According to Scheps et al. [7] their data at long wavelengths are noisy. The long 
range interactions should be reliable [17] and we have explored the effects of modifying the X 2 S potential in the short 
range repulsive region. However we were unable to obtain precise agreement with the measurements. 

In Fig. 9, we present the red and blue wings of the line profiles in the high pressure limit for temperatures between 
200 K and 3000 K. The total energy emitted per second in the red wing at wavelengths longer than 690.8 nm, given 

by 

ffo-Ai 

h = \ hvk v dv, (12) 

Jo 

where Ai corresponds to a wavelength shift of 20 nm red from the line center, and the total energy emitted per second 
in the blue wing at wavelengths shorter than 650.8 nm, given by 

/>oo 

I2 = / hvk l/ dv 1 (13) 
J i/ +A 2 

where A 2 corresponds to a wavelength shift of 20 nm blue from the line center, are listed in Table IV for temperatures 
between 200 K and 3000 K. 

The absorption coefficients are shown in Fig. 10. They arc similar in shape to the semi-classical calculations of 
Jungen and Staemmler [8] and of Erdman et al. [9]. There is a satellite near 536 nm which also appears in the 
emission spectrum but otherwise the absorption coefficients decrease steadily with separation from the position of 
the resonance line. There is an indication of a satellite near 600 nm in the calculation of Erdman et al. [9] based on 
different potential curves which may have the same origin. 

VI. CONCLUSIONS 

We have carried out quantum mechanical calculations for the Li-He emission and absorption spectra at temperatures 
from 200 K to 3000 K and wavelengths from 500 nm to 1000 nm. We find a blue satellite near 536 nm for T = 3000 K 
and a red satellite near 870 nm for T w 200 300 K in the emission spectra, and a blue satellite near 536 nm 
for T w 2000-3000 K in the absorption spectra. At 670 K, our emission coefficients are in good agreement with 
experiment. 
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TABLE I: Vibrational energy levels of the A 2 II state of the LiHe molecule. 



v 


1 


2 


3 


4 


5 


6 


E v,j=u ( cm -i) _837.7 


-549.3 


-327.7 


-168.2 


-67.11 


-15.69 


-0.9170 



TABLE II: Resonance energies E^.' J and widths r f r ' J of the quasi-bound levels of the A 2 IT state. 



J 


V 


7—171 / /TT j_ \ 

Ep J (Hartree) 


7^ ,J (Hartree) 


2 


1 


2.14x10"' 


2.1xl0~ M 


6 


1 


1.34xl0~ 5 


6.5xl0~ 7 


9 


1 


1.04xl0~ 5 


7.1xl0~ 13 


10 


1 


6.56xl0~ 5 


1.2xl0~ 6 


11 


1 


1.19xl0~ 4 


1.8xl0~ 5 


13 


1 


9.15xl0~ 5 


2.3xl0~ 9 


14 


1 


1.94xl0~ 4 


2.1xl0~ 6 


15 


1 


2.91xl0~ 4 


2.4xl0~ 5 


16 


1 


1.23xl0~ 4 


1.3x10 12 


17 


1 


2.90xl0~ 4 


7.0xl0~ 8 


18 


1 


4.50xl0~ 4 


6.0xl0~ 6 


19 


1 


1.82xl0~ 4 


3.6xl0~ 15 


20 


1 


4.21xl0~ 4 


1.6xl0~ 9 


21 


1 


6.57xl0~ 4 


7.0xl0~ 7 


22 


1 


2.84xl0~ 4 


6.2xl0~ 17 


22 


2 


8.81xl0~ 4 


1.7xl0~ 5 


23 


1 


6.07xl0~ 4 


3.4X10" 11 


24 


1 


9.31xl0~ 4 


4.2xl0~ 8 


25 


1 


1.25xl0" 3 


3.0xl0~ 6 


26 


1 


1.55xl0" 3 


3.1xl0~ 5 
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TABLE III: Lifetimes (ns) of selected ro-vibrational levels (v, J) of the A 2 n state. 





w = 


v = 1 


v = 2 


v = 3 


v = 4 


v = 5 


J = 


59.2 


51.5 


44.1 


37.7 


32.4 


28.6 


J = 5 


58.8 


50.9 


43.5 


37.0 


31.7 


27.7 


J= 10 


57.4 


49.4 


41.8 


35.0 






J = 18 


53.0 


43.7 
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TABLE IV: Total energy emitted per second (ergs cm 3 s 1 ) denned in Eqs. (12) and (13) for the red and blue wings, 
respectively. 



T(K) 


Red wing 


Blue wing 


200 


4.68x10"^ 


1.46x10" 


2y 


300 


7.59xl0~ 26 


9.79x10" 


29 


670 


1.64xl0~ 26 


8.27x10" 


28 


1000 


1.22xl0~ 26 


1.59x10" 


27 


2000 


9.39xl0~ 27 


3.51x10" 


-27 


3000 


8.67xl0~ 27 


4.88x10" 


-27 
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FIG. 1: (a) Adopted potentials V(R) in atomic units for the X 2 S, A 2 Il and B 2 E states, (b) Difference potentials V A 2 n (R) — 
V X 2- E {R) and y B 2 E (J?) — V X 2 S {R) in atomic units. 
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FIG. 4: Spontaneous transition probabilities W v , Eq. (11), in atomic units for vibrational levels v = — 5forJ = (upper 
panel), 5 (upper middle panel), 10 (lower middle panel) and 18 (lower panel). 
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FIG. 6: Spontaneous transition probabilities W v , Eq. (11), in atomic units for ro-vibrational levels (v, J) = (5,0) and (5,5). 
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FIG. 8: Comparisons of our calculated emission coefficients with experimental measurement of Scheps et al. [7] and previous 
theoretical calculations of Herman and Sando [10]. For the blue wing, the high and low pressure curves are identical. Unit gas 
densities are used, riLi = njj e = 1 cm" 3 . 
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